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Abstract 


Strong deformation of ultrashort laser pulse shapes is unavoidable when delivering high inten¬ 
sities at remote distances due to non-linear effects taking place while propagating. Relying on 
the reversibility of laser filamentation, we propose to explicitly design laser pulse shapes so that 
propagation serves as a non-linear field synthesizer at a remote target location. Such an approach 
allows, for instance, coherent control of molecules at a remote distance, in the context of standoff 
detection of pathogens or explosives. 
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I. INTRODUCTION 


Manipulating atoms and molecules using quantum control or pump-probe excitation 
schemes requires precisely tailored electromagnetic fields jl-5] which require stability in 
time and space. Most of the manipulations rely on non-linear interactions of broadband 
pulses with molecules, and therefore high peak intensity (GWcm -2 to tens of TWcm -2 ) 
and ultrashort (femtosecond) pulses. Therefore, the pulse propagation from the laser source 
to a remote target location is likely to encounter non-linear spatio-temporal perturbations 
as, e.g., Kerr-induced wavefront distortion, self-phase modulation and temporal pulse shape 
alteration, in addition to linear perturbations such as diffraction, dispersion and aberrations 
that stem from air turbulence. In laboratory-scale experiments, these perturbations can be 
handled using deformable mirrors and temporal or spectral pulse shapers. Linear deforma¬ 
tions can be corrected over large distances as well. This is routinely performed with adaptive 
wavefront corrections based on artificial guide stars for large astronomical telescopes [6], or 


0 . 


the correction of turbulence perturbations with temporal pulse shapers 131 . These corrections 
can be either static or dynamic, e.g., based on a feedback loop ji|. 

In contrast, non-linear spatio-temporal phase shifts accumulate and give rise to dramatic 
modifications of the laser waveforms. If Kerr self-focusing overcomes diffraction, which re¬ 
quires the peak power to exceed a critical power of 3 GW, the beam can even experience 
a catastrophic collapse. The arrest of the_collapse by self-defocusin g pr ocesses (ionization 


or higher-order optical non-linearities 


BQ) 


leads to filamentation 


121-15]. The complex 


propagation dynamics associated to this highly non-linear propagation mode results in a 
dramatic reshaping of the pulse, challenging any explicit design of the initial pulse required 
for open loop approaches. Furthermore the lack of feedback from a remote pulse character¬ 
ization at the target location prevents a closed-loop approach. 

Active remote sensing, including standoff detection of explosives or harmful molecular 
agents, greatly benefits from non-linear interactions induced by ultrashort high intensity 
lasers. In particular, coherent Raman spectroscopy 4 Il6l~l22|. multi-photon excited fluo¬ 


rescence 


23H25], and even atmospheric lasing [26 29] require a precise and stable spatio- 


temporal waveform at a defined distance. A particularly perturbation-sensitive strategy 
is the discrimination of almost identical molecules using coherently controlled multiphoton 
30-24], like optimal dynamic discrimination (ODD). In order to transpose 


fluorescence 
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the successful demonstration of bacteria or biomolecules discrimination from the labora¬ 
tory to an atmospheric standoff detection workframe, non-linear phase perturbations have 
to be included when manipulating optical waveforms: If the B-integral (B = k 0 f* n 2 Idz 
12| ) characterizing the accumulated non-linear phase shift exceeds a fraction of it, these 
perturbations have to be pre-compensated. Since the required laser power often exceeds 
the critical power in order to propagate over large distances, hlamentation is also likely to 
occur. The phase pre-correction therefore requires a precise knowledge of this self-induced 
spatio-temporal phase modulation in filaments. 

In the present work, we propose to make use of the unique properties of hlamentation 
(self-broadening of the laser spectrum, self-cleaning of the spatial mode, temporal effects like 
self steepening,..) jl2-l^] as part of the waveform synthesizer required to control molecules 
at a target distance, instead of trying to prevent it. Indeed, hlamentation allows delivering 
extremely high intensities remotely, circumventing the diffraction limit 0 . It is therefore 
well suited to initiate, control, and optimize non-linear processes, that require high-intensity, 
shaped pulses on the target. 

Such approach requires to tailor a pulse so that it acquires the desired shape at the remote 
location of interest, i.e., during or after hlamentation. In spite of intensity clamping H 
and of the extreme phase sensitivity at the non-linear focus j^Z], we recently showed that 
hlamentation is reversible jsS]. It is then possible to numerically back-propagate a target 
electric held to the initial conditions giving rise to it. Here, we show that this previously 
unexpected reversibility can be applied to explicitly design initial pulse shapes which yield 
the required target electric held for specihc applications after hlamentation or, more gen¬ 
erally, after non-linear propagation. We illustrate this approach and its versatility by two 
examples featuring the shaping of temporal intensity, and spectral phase, respectively. 


II. NUMERICAL METHOD 


A. Non-linear propagation of ultrashort laser pulses 


The reversibility of laser hlamentation relies on the inversion of the propagation direc¬ 
tion in the governing equations[38|. i.e., the non-linear Schrodinger equation (NLSE), or 
the unidirectional pulse propagation equation (UPPE)0|. These equations describe the 
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propagation of the complex electric field £ along the axis z (assuming axisymmetry): 

d z e = \[k z - )e + — 

v g c z k ; 

where v g is the group velocity, e and m e the charge and mass of the electron, respectively, 
k z = y/k 2 (a;) — k^_, with k(u) being the wave vector and k± its transverse component, and 
n 2 the Kerr index. t(lu) = (u en + i u)/(v 2 n + u> 2 ), u en is the collision frequency between free 
electrons and neutral atoms and the tilde denotes a Fourier-Hankel transform: 
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( 1 ) 


£(k±,z,co ) = / / rJ 0 (k±r)£(r, z,t)e lut dtdr, 


( 2 ) 


where Jq is the zeroth order Bessel function. The free-electron density p evolves as: 


d t p = W(\£\ 2 )(N-p) + ^\£\ 2 -g(p), 


( 3 ) 


where W(|£| 2 ) describes the probability of ionization calculated with the Perelomov, Popov, 
Terentev (PPT) formula jd^], a is the inverse bremsstrahlung cross-section, U\ the ionization 
potential, N is the number density of molecules, and g is the recombination function. The 
last term in equation (??) accounts for ionization-induced losses, and is calculated as 


m- U ^N- p) e. 


( 4 ) 


At 400 nm, plasma constants are given in ref. 
800 nm, these constants are given in ref. 


12] and the Kerr index in ref. |4lfl. At 


B. Back-propagation of filamenting pulses 


The pulses bearing the desired shapes (as determined from prior theoretical arguments 
or laboratory experiments) at the target location are back-propagated by changing dz to 
—dz in the propagation equation (??), as detailed earlier [28]. 

Reversibility however neither ensures that any arbitrary pulse shape can readily be back- 
propagated, nor that it will follow a filamentary path during its back-propagation. In other 
words, designing an initial pulse shape that will synthesize a target waveform after under¬ 
going filamentation requires specific features. For instance, to back-propagating arbitrary 
pulses can lead to intensity divergence due to the gain induced by the inversion of the loss 
term related to ionization 0 - 
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FIG. 1. (Color online) Principle of pulse hybridization for designing remote, post-filamentation 
pulse shapes for standoff detection applications, (a) Local a priori calibration experiment defining 
the target pulse shape; (b) Forward filamentation propagation defining a pulse framework able to 
yield filamentation; (c) Back-propagation of the hybridized pulse, that mixes the two previous ones 
and defines the required initial pulse shape; and (d) Remote optimal experiment. The forward 
propagation of the hybridized pulse yields a shape identical to the target on the beam axis, thanks 
to the numerical stability of the filamentation reversibility. 


C. Pulse hybridization 

Filamentation is likely to occur for high intensities and long propagation distances so that 
strong coupling between the spatial and temporal phase takes place. For this reason, we had 
to specifically design hybrid pulses encompassing the drastic spatial deformation induced by 
filamentation. To design such pulses, we developed a hybridization strategy. As sketched in 
Figure [lj this is achieved by embedding the desired target shape (Panel a) into the center of a 
post-filamentation propagated pulse (Panel b), therefore resulting in backward filamentation 
(Panel c) while bearing the tailored shape. The size of the transition region has to be defined 
as a trade-off between three constraints: (i) ensure that the central region of the pulse is 
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adequately shaped, (ii) provide a filamenting reference over a peripherial region bearing a 
substantial part of the pulse energy, and (iii) offer a sufficient intermediate zone to ensure a 
smooth transition. 


More specifically, hybridization was performed in two distinct ways for temporal intensity 
and spectral phase shaping, respectively. In the former case, we applied the target electric 
field e t (Figure dh) at the center of the beam profile, and progressively mixed it over the 
central region (r < r m = 0.75 mm) with a reference field (e p of phase <f> p ) originating from 
the forward filament propagation (Figure HJ)): 




where the sine weighting was chosen so as to ensure a continuously differentiable radial tran¬ 
sition between the two pulse components. Beyond r m , the propagated pulse amplitude e p 
and phase <h p are used to maintain the broad spectrum generated by the forward-propagated 
filament. Such hybridization ensures the stable back-propagation as a filament. Simultane¬ 
ously, the post-filamentation shape is preserved in the outer region of the profile. It accounts 
for 82% of the pulse energy, ensuring the occurrence of filamentation in the backward prop¬ 
agation. 

In the case of spectral phase shaping, we aimed at preserving simultaneoulsy the ra¬ 
dial continuity of the electric field envelope in the temporal domain, and the wavefront 
curvature at each wavelength and for each radial position. We therefore hybridized the 
phase rather than the complex electric field. More precisely, we considered a propagated 
(post-filamentation) pulse of phase <fi p (r,uj). At each wavelength, the difference between the 
propagated phase on the side (0 p (r, u>)) and the center (</> p (0,u;)) describes the radial phase 
dependence, and the associated wavefront curvature. The desired phase <fi t (u;) was added to 
this curvature, defining in the Fourier space the target pulse to be backpropagated as: 



As hybridization affects the peak power and energy of the pulse, the latter was adjusted so 
as to recover a power allowing backward filamentation without triggering intensity divergence 
by the backward gain. 


The need to hybridize the electromagnetic waveform is illustrated in Figure Eh. An 


800 nm, 0.44 mJ non-hybridized shaped pulse, with a purely spherical wavefront curvature 
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(blue dash-dotted line) does not enter a filamentation regime when back-propagated. In 
spite of a peak intensity of 70 TW cm -2 in the focal region, self-guiding does not occur, as 
evidenced by the short length of both the high-intensity region (5 cm) and the ionized region 
(6 cm above 10 14 cm -3 , see inset), similar to the double Rayleigh length (5 cm) associated 
with a Gaussian pulse with identical beam diameter and focusing. A similar, filament-free 
backward-propagation dynamics is also obtained with longer focal lengths (/ = 2 and 3 m). 
This behavior contrasts in the case of back-propagation of a properly hybridized pulse (black 
line), the intensity of which is clamped at the same intensity of 70 TW cm -2 , showing that 
it is guided over 20 cm, and ionized over 17 cm. 


III. RESULTS 


A. Remote synthesis of pulses with temporal intensity control 


We first investigated the design of a pulse that will provide, after propagation, a desired 
temporal intensity shaping. We simulated the remote generation of a train of three Gaussian 
sub-pulses of 12 fs duration (FWHM) at 800 nm, with 32 fs separation.This kind of structure 
is typical, for instance, for stimulating vibrational Raman scattering jd^]. This target pulse 
was back-propagated to z = 0 in order to define the initial pulseshape required to generate 
this target triple pulse after filamentation. 

The target and back-propagated input pulses are notably different (Fig. [BJ). The spec¬ 
trum shrinks during the back-propagation, narrowing down to the spectral support for a 
Fourier Transformed 30 fs pulse envelope, while the temporal shape simplifies into a single 


peak envelope with a shoulder. Accordingly, the spectra! 
synthesized experimentally by a standard pulse shaper 


phase is sufficiently smooth to be 


431 ]. The manipulation can there¬ 


fore be seen as a controlled pulse splitting to generate a train of 3 Gaussian subpulses at 
the target location. The error between the target pulse (Fig. EK) and the pulse obtained 
by numerically forward-propagating again the back-propagated pulse (Fig. [3b) stays below 
0.2 %, suggesting that the proposed non-linear pulse manipulation could be implemented 
experimentally. 

The significant reshaping of the pulse is related to the complex dynamics occurring dur¬ 
ing filamentation (Figure Eb>): Between z = 0.85 and 1.05 m, the intensity reaches up to 
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FIG. 2. (Color online) Back-propagation of a target triple pulse, (a) On-axis intensity and electron 
density ; (b) Evolution of the temporal dynamics of the pulse along the propagation distance. 

70 TWcm -2 , the pulse maximum shifts temporally, and then splits into three sub-pulses, 
two of them developing in the far-held during the post-hlamentation propagation stage. This 
complex evolution of the pulse illustrates the rich dynamics of hlamentation, and the asso¬ 
ciated difficulty of explicit pulse design. It also demonstrates the relevance of the proposed 
pulse hybridization procedure: It adequately prepares pulseshapes such that hlamentation 
acts as a complementary nonlinear pulse shaping step. From a more fundamental point of 
view, in this example, tailoring only 18% of the pulse energy has a crucial inhuence on the 
hnal pulse shape, as expected due to the high sensitivity of non-linear propagation to the 
initial conditions. 


B. Remote synthesis of pulses with spectral phase control 

The temporal hybridization used in the above example, although intuitive, leads to an 
inhomogeneous temporal prohle along the radius of the back-propagated pulse (Figure [3b), a 
shape that cannot be realized experimentally with today’s technology. In contrast, the phase 
hybridization can be expected to result in a more homogeneous temporal shape along the 
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FIG. 3. (Color online) Tailoring a triple pulse at 800 nm after laser filamentation. (a) Hybridized 
target pulse, featuring a triple pulse around its center, surrounded by the output of a regular 
filamenting pulse; (b) Corresponding on-axis spectrum and spectral phase on the beam axis; (c,d) 
Same data for the back-propagated pulse, i.e., the pulse shape that has to be propagated to get the 
target shape. 


radius. Furthermore, phase shaping is the most widely used approach for coherent control, 
e.g. in view of optical discrimination of molecules featuring similar spectroscopic properties 

With highly selective standoff detection applications in mind, we investigate the possibil¬ 
ity of discriminating the presence of molecules relevant to bioaerosols: Riboflavin (RBF) and 
flavin mononucleotide (FMN) ( 22 ]. This demonstration was successfully performed in the 
laboratory { 3 I, 33], in which the optimal waveforms were determined. Figure UK illustrates 
the complexity of this optimal pulse, featuring seven significant sub-pulses with different 
chirps and temporal shapes. Realizing this biomolecular identification at a distance implies 
delivering these high-intensity pulses at 400 nm with the spectral phase maximizing the 
depletion of fluorescence from RBF as compared to that of FMN (and vice-versa). 

The contrast between target and back-propagated pulses is clearly observed on both the 
spectrum and the spectral phase. The target and initial pulses exhibit differences up to 60 
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FIG. 4. (Color online) Tailoring of an optimally phase-shaped pulse for remote selective detection 
of RBF with regard to FMN: Wigner plots (with 12 fs FWHM Gaussian gate) of (a) the target 
pulse shape, and the central part of (b) the target and (c) the back-propagated pulses; (d) Radial 
distribution of the intensity of the back-propagated pulse. 


% for the central peak intensity, and up to 0.9 ir for the dephasing of the spectral phase 
(Figure Eb,b). In the temporal domain (Figure Eb), this results in a splitting of the central 
peak of the back-propagated pulse (between t = -100 and 60 fs). These deformations are 
also visible in both the temporal and frequency spaces, as evidenced by the Wigner plot 
at the center of the pulse (Figure |Ub,c). While the target pulse has a double structure 
without chirp (Figure 0b), the above-mentioned splitting results in three sub-pulses in the 
back-propagated pulse (Figure 0b). Furthermore, the two first ones have opposite chirps. In 
the temporal domain, the pulse reshaping mainly occurs on phase rather than on intensity, 
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FIG. 5. (Color online) Target and back-propagated waveforms on-axis (a) spectral phase, (b) 
spectrum and (c) temporal pulse shapes of an optimally phase-shaped pulse for remote selective 
detection of FMN with regard to RBF. 

as evidenced by Figure Eh. 

The shape of the back-propagated pulse depends much less on the radial position than 
with the temporal hybridization (Fig Uhl). This homogeneity appears even more clearly when 
considering that the beam is more than twice as large as in the previous case, and that the 
phase curvature is very small on the first two millimetres around the axis. 

The spectacular reshaping of the pulse during the back-propagation is again related to 
the phase scrambling induced by non-linear propagation associated with filamentation. As 
illustrated on Figure IHl the back-propagation of our target pulse displays filamentation 
over almost 10 cm, as evidenced by the ionized section and the refocusing cycle. This is 
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FIG. 6. (Color online) (a) Intensity and (b) free electron density dynamics during the back- 
propagation of a spherically focused pulse at 400 nm, and of a hybridized filamenting pulse op¬ 
timizing the discrimination of RBF vs. FMN. 

again made possible by preparing an adequate pulse hybridization. We compared these 
results with the back-propagation of a 0.84 mJ pulse with the same spectral phase but a 
spherical wavefront curvature instead of pulse hybridization ensuring filamentation during 
the back-propagation. The non-hybridized pulse does not enter a filamentation regime when 
back-propagated, although it reaches the same peak intensity of 76 TW cm -2 in the focal 
region (1.6 m, Figure E|). 


IV. DISCUSSION 


The results unambiguously show the possibility of explicitly tailoring the output of laser 
filaments in air to a prio ri arbitrary waveforms. Potential applications range from selective 


standoff detection 


0,ll6 


16d21 


23h27| to lightning control [44], atmospheric water condensation 
|4. r )l | . and radiative forcing modulation [4cl | . Similar techniques could in principle be applied 
to any other remote application. It only requires that the target shape can be explicitly 
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designed and/or empirically determined a priori , e.g., via a closed-loop algorithm or other 
optimization techniques in laboratory (local) experiments beforehand, or from numerical 
simulations of the phenomenon of interest. For example, we successfully back-propagated 
several other target waveforms, including Gaussian pulses, double pulses, or phase controlled 
pulses maximizing FMN vs. RBF fluorescence depletion instead of minimizing it. 

The experimental implementation of the pulse tailoring technique proposed in the present 
work requires the back-propagated pulse shapes to be compatible with present technical 
constraints, especially a decoupling of the spatial and spectral phases. Our work shows that 
the way to perform the pulse hybridization is critical in that regard. 

Finally, the proposed pulse design technique is based on the reversibility of filamentation, 
that does not depend on the specific filamentation or ionization model |38]. Consequently, 
our technique is also independent on these models. 


V. CONCLUSION 

In summary, we have demonstrated the explicit remote tailoring of arbitrary pulses after 
non-linear propagation, and more specifically at the end of laser filamentation. It opens 
the way to applications requiring the remote delivery of specific laser waveforms with high 
intensity, in order to optimize atmospheric applications like lightning control or laser-induced 
condensation, or selective remote sensing. 
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